Quantum Non-Equilibrium Steady States Induced by Repeated Interactions 
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We study the steady state of a finite XX chain coupled at its boundaries to quantum reservoirs 
made of free spins that interact one after the other with the chain. The two-point correlations 
are calculated exactly and it is shown that the steady state is completely characterized by the 
magnetization profile and the associated current. Except at the boundary sites, the magnetization 
is given by the average of the reservoirs' magnetizations. The steady state current, proportional 
to the difference in the reservoirs' magnetizations, shows a non-monotonous behavior with respect 
to the system-reservoir coupling strength, with an optimal current state for a finite value of the 
coupling. Moreover, we show that the steady state can be described by a generalized Gibbs state. 



Understanding non-equilibrium behavior of quantum 
systems on the basis of general principles is one of the 
more challenging prospects of statistical physics. In par- 
ticular, the complete characterization of the so-called 
Quantum Non-Equilibrium Steady-States (QNESS), i. e. 
stationary current full quantum states, is of primer and 
central focus since they are possible candidates for play- 
ing a role similar to Gibbs states in constructing a non- 
equilibrium statistical mechanics [I, 2, 3, 0, 0, 0, S Si • 

To elucidate the general guiding principles for a non- 
equilibrium statistical mechanics, exactly solvable mod- 
els play a central role. Among many models, the XX 
quantum chain is one of the simplest non-trivial many- 
body system. Its A'^-sites Hamiltonian is given by 



Hs = 



2 ^ 

1=1 



N 

1=1 



(1) 



where the as are the usual Pauli matrices, J is the 
exchange coupling and h a transverse (possibly exter- 
nal) magnetic field. Due to the fact that its dynamics 
can be described in an explicit way, the one-dimensional 
XX model has been extensively studied in various non- 
equilibrium contexts 0. On the experimental side, the 
most promising perspectives for these studies come from 
the ultra-cold atoms community since the XX model ^ 
can be mapped on a one-dimcnsional Hard-Core boson 
(Tonks-Girardeau [13]) model through the transforma- 
tion = (erf -I- icrf)/2 and 5/ — {af — iaf)/2. Experi- 
ments on such ID hard-core bosons have been performed 
with Rubidium atoms in both continuum [11] and lattice 
[l^ versions. 

Antal et al. [l^] studied the ground state of the Ising 
and XX chain Hamiltonian with the addition of a mag- 
netization or energy current via a Lagrange multiplier. 
The ground state of the effective Hamiltonian Hs — \ J 
was interpreted as a non-equilibrium stationary current 
full state. Such an effective Hamiltonian was supposed to 
capture locally the essential features of a finite chain cou- 
pled at its boundary sites to quantum reservoirs. Soon 
after, Ogata [3], Aschbacher and Pillet [is'] considered 



the anisotropic XY steady state induced by the unitary 
dynamics, V{t) — e~'*^^, starting with an initial state 
in which the left and right halves are set at inverse tem- 
perature /3l and They showed that the QNESS 
can be effectively described by a generalized Gibbs state 
~ ^-0Hs+SY ■^yjjej.g ^ = {fiL+fin) /2 is the average inverse 
temperature, 5 = {(3l — Pr)I'^ is a driving force coupled 
to a long-range operator Y commuting with Hs ■ The op- 
erator Y is given by y = X^z^i A*'^ where Yi are currents 
operators associated to ^th sites conserved quantities and 
where the coefficients /i/ show a power law decay (~ 1/0- 
The long range of F is a signature of a strong nonlocal 
properties of the QNESS which is believed to be a generic 
feature of NESS [ij] . One may notice that the Antal et 
al. steady state relates to the same effective Hamiltonian 
truncating the current series Y — ^{Yi to the first few 
terms. 

In this work we study the A^-site isotropic XX-chain 
(H]) coupled at its boundary sites to quantum reservoirs at 
different temperatures. In the anisotropic case a recent 
study has been reported in [3] with Markovian baths. 
Here, the left and right reservoirs are made of an infi- 
nite set of non- interacting spins 1/2 in the same trans- 



verse field h with Hamiltonian H 
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The system-reservoir couplings are 



implemented via a repeated interaction scheme, meaning 
that each subsystem (particle) composing the reservoirs 



are interacting with the system one after the other 17]. 
To have a physical picture of such an interaction scheme, 
one may think of a laser (or particle) beam falling on the 
system. The system-reservoir interaction is given by the 
time-dependent Hamiltonian V{t) — VL(t)-f yR(t) where 



L{R) 



el 



l)r, tit] with 



Vl^r) = - + Ki.N)bn^"-^] ' (2) 

selecting the nth left and right reservoir spins in the time- 
interval t €](n — l)r, nr]. 

We start at t = with a system-environment decou- 
pled initial state p(0) — ps{0) ® Pe{0), where ps is 



2 



0" ^^ 



FIG. 1: Sketch of a system in contact at both ends with 
repeated interaction reservoirs. 



an equilibrium state of the system and where the en- 
vironment density matrix is given by 
(^N'iPn ^ Pn) ^ith one-particle thermal density matri- 



ces Pn 



1 — m 



L(R) 



-)(-| with cr^|±) = ±|±) [18|. At t = 0+, the 
first left and right reservoir spins start to interact with 
the left and right system boundary spins for a time r 
through the hopping term Vl^j^y At f = t+, the first 
reservoir spins are replaced by the second ones interact- 
ing with the system through for a time r. The 
process is then repeated again and again [l^. Iterating 
the process n-times, the reduced density matrix associ- 
ated to the system part is given at time t = nr hy 

psinr) = TrE{e-'^"T':lp((n~l)T)e''-"T"ol^ 

= rr„(c/j"Vs((n-l)r)®p„[/j"^H (3) 



where H, 



{"} 



^Tot - Hs + He 

in the time-interval [(n— 1)t, nr] with the interaction part 
y(n) ^ Y(n)_^y(n)^ and whcrc Uj"^ = (ffs+v/">+ff„) 
with Hn = H^+ H^. 

To solve the recursive dynamical equation we in- 
troduce the fermionic representation of the coupled part 
of the total Hamiltonian, i7^^"> = Hs + v/"^ + H„, 
which is of the form ([T]) with N + 2 sites, via the usual 
Jordan- Wigner transformation |20|: 



V^'^^ is the total Hamiltonian 



where the Fs are Majorana real (Clifford) operators sat- 
isfying Ft = F and {Ff,F^} = 25ij6a(3. Notice here that 
the zeroth label is associated to the left-reservoir spin 
and the N + 1th to the right one, keeping the l,..,iV 
labels for the system. The interacting part of the to- 
tal Hamiltonian takes the form H^^'^'^ = l/4rtTs„r 
where pt = (F^, F}, F^^+i, T^, F^+J is a 2{N + 2)- 
component operator and where the 2{N -I- 2) x 2(A^ + 2) 

matrix T is given by Tsn — ( ) '^ith the tridi- 

V ^Sn / 

agonal matrix (Cs„);„ = -i{h5im + JArn-i + JnAm+i) 
containing system-environment couplings Jq — Jn = Je 



and J; = J VZ = 1, ...,N — 1 for the system part. The 
interesting point in using the Fs is that their time evolu- 
tion, generated by H^^^'^\ is simply given by a rotation 
R{t): T{t) = Rit)T(0) = e-^*^S"r(0) li. 

Since the total initial state is Gaussian in terms of 
fermions, the reduced system density matrix remains 
Gaussian during the repeated- interaction process 21 , 2^ ■ 
No many-body interactions are generated during the time 
evolution. As a consequence, thanks to Wick's theorem, 
one may characterize completely the state of the system 
at any time by its two-point correlation matrix Gs{t) 
defined by {Gs{t))^, = iTrs {[iTs)k. {Tshjpsit)} . 
Along the same lines, the two-point correlation matrix 
G'f'^"''(nr) characterizes completely at time nr the total 
state of the nth environment copy -I- system. For these 
reasons, instead of computing directly the system-density 
matrix, thanks to Q, we will study the correlation ma- 
trix and reconstruct afterwards from it the density ma- 
trix. 

Ordering the Ft = (r^,rg) such that the first part, 
= (r|^, r^), is associated to the nth copy of the en- 
vironment and that contains the components of the 
system, we write the rotation matrix R(t) as 



Re Res 
RsE Rs 



(5) 



Using this decomposition into ^ , one arrives at the fun- 
damental dynamical equation governing the system: 

Gsinr) = RsGsiin - 1)t)Rs^ + RseGeRse^ ■ (6) 

The 4x4 environment correlation matrix Ge describes 
the initial environment state. It is evaluated with re- 
spect to the environment two-spin initial state Pn ^ Pn- 
If the system is decoupled from the reservoirs {Je = 0) 
the rotation matrix splits into a block diagonal form, 
with Res = 0, Rse = 0, reflecting the separate unitary 
evolution of the system and the environment through 
Rs and Re respectively. Introducing the infinitesi- 
mal generator C through the dynamical map agi^X) = 
e~'^^{X) = RsXRs^ one may iterate equation ^ which 



(^) becomes in the continuum limit Gs{t) 



-tc 



(Gs(0))- 



J^ds e-'^ (Ge^ with Ge = RseGeRIe- The explicit 
form of the generator C depends on the proper way 
one rescales the interaction couplings in the continuum 
limit T —^ 0. To take into account the non-trivial ef- 
fect of the interaction between the system and the en- 
vironment, one has to rescale the interaction couplings 
Je — ^ Je/^/t- Other rescalings give either trivial lim- 
its or no limit at all 17, 2^. The total T- matrix takes 

Te e/V^ \ 
et/v7 Ts J 



then the form T; 



Sn 



where Te 



and Ts are T-matrices of the environment and system 
parts respectively while O is a 4 x 2N matrix describ- 
ing the system-environment interaction with components 



3 



(6)fei = -iJE{Ski5iN+i - Sk2Sii + SksSnN - SkiSm)- De- 
veloping the exponential R{t) = e~''^"^^" to the lowest 
order in r, and projecting to the system part leads to the 
Linblad- like differential equation i9tGs(f) = ~C{Gs{t)) + 
etG^e with the generator £(.) = i[Ts, .] + ^{Q^Q, ■}■ 

Using the expressions of the interaction matrix Tsn 
and the antisymmetry of the correlation matrix, one ar- 
rives finally at two coupled N x N matrix equations: 



dtGo = i[Gd, Cs] - #Ao - JIMe 



(7) 



where the NxN matrices Gd and Go are defined through 

Gs = ( ^° I sind where Cs is the restriction of 

Gsn to the system part. The matrix Me with elements 
{ME)ki = TnL5ki5ii mjiSkN^m contains the left and 
right environment magnetizations. The relaxation ma- 
trices ^-d,o^ related to {6^9, Gs}, are given by 



/ 



A. 



G 



V 



Ctv-i 1 
'2'Gcni 






Gcn2 



^Ci jv-i 






\ 



(8) 

with C = 0,0?. Notice that these matrices contain an 
identically vanishing internal (TV — 2) x {N — 2) square. 

It can be proven that the steady state is unique [3, 
and reached exponentially with a relaxation time T ^ 
[3]. The steady state correlation matrix G*g obeys 
([7]) with the left-hand side set to zero. In the vanish- 
ing square sector of Ad,o, thanks to the antisymmetry 
of the correlation matrix Gs, one derives from ([7]) the 
space translation invariance of the matrix GJ: {G*j)ki — 
G*^{1 — k). It implies in particular that the steady state 
magnetization current, {J^)* = 2JGl,^ = 2Jj*, de- 
fined through the quantum continuity equation cr| — 
i[Hs,crk] = J^k-i ~ J^™' is constant all along the chain. 
Using this translation invariance, one reduces the full set 
of steady equations to 



rriL 



= 73 

f hi ■ 



-- m 
m' 



N 



(9) 



with 7 = 2J ^^'^ where — (a^) — —{G*)kk is the 
steady state magnetization at site k, which takes a con- 
stant value denoted m* Vfc — 2, ...,7V — 1. One finds 
that all other correlation matrix elements are identically 
vanishing in the steady state. From ^ it appears that 
the four unknowns j* , m* , m| and are functions of 7 
and niL.R (that is on /3l,_r and h) only and consequently 
size- independent. Solving © one finally finds that the 
exact steady state properties are fully characterized by a 
stationary current 



7 rriL-mR 




FIG. 2: Rescaled stationary current as a function of the cou- 
pling ratio 7 = iS-. 



and a flat magnetization profile 



TO,. 



rriL + mR 



V/c 



,N-l 



(11) 



for the bulk spins and m\ = to* -I- 7j* , = m* — 7j* 
for the boundary sites. The size-independence of these 
quantities has to be related to the perfect ballistic nature 
of the elementary excitations transport properties. One 
may notice that the boundary values are deviating from 
the flat profile by an amount which is proportional to the 
current value. However, in the large reservoir coupling 
limit, 7 — > 00, the magnetizations of the left and right 
boundary sites tend to the corresponding reservoir val- 
ues. It is interesting to note that while the bulk magne- 
tization profile is independent of the interaction strength 

ratio, 7 = the current show a non-monotonous be- 
havior, see figure with a maximal current state for 
7 = 1. One may explain this behavior by noticing that 
at small 7 the system is very weakly coupled to the reser- 
voirs and it is very unlikely to inject a particle (or flip the 
boundary spin) at the boundary site, leading to a small 
current value. On the contrary for large 7 the coupling 
to the reservoirs is much larger than the chain coupling 
and it is very easy to flip the boundary spin, but hard to 
propagate this flip along the chain and this leads again 
to a small current. 

From the knowledge of the two-point correlation ma- 
trix one can deduce the steady state of the system, 
which will appear to be of generalized Gibbs form. To 
show that, let us introduce Dirac Fermions: F^ = + Cfc 
and F^, = i{c^ — Ck) and search for a quadratic form 
Q = ^ • j cl AijCj such that p*g ^ e~^. The coupling ma- 
trix A, giving Q, is deduced from the correlation matrices 
G^ o '-'^ equivaientiy, using the Dirac Fermions represen- 
tation, from Afij = Tr{c^Cip*g} thanks to the relation 



^ l-7V ^ 



1 + 7' 



A = ln 

In the large system size limit, after the diagonalization 
of the correlation matrix A/", one obtains from the pre- 
vious relation Ak+i,k = (-l)'^fc,fc+i = where ai^o = 



(10) sgn'(/)^ \(^iyzi(Ill)- zi(k:3l)] with z{x) = 



4 



- 1 and = In [z{^) / z{^^yj . Using that into 
Q one finally obtains the generalized Gibbs form p*g ~ 

g-aoQo/2-E,^o".Q< with = C++,Cj + (-l)'c+c,+, a 

set of conserved quantities. 

To interpret this result, and eventually extract an ef- 
fective temperature characterizing the QNESS, consider 
first the undriven situation j* = for which I3l,r — P 
(^L,fl = «*)• In that case, since ai oc {j*f, all a/^o = 
while ao reduces to In = f3h. Consequently the sta- 
tionary state reached in our setup is described by the 
equilibrium Gibbs state e~^^ with a free spin reference 
Hamiltonian = ^AI^ = ^Qo + const, and a temper- 
ature f3~^ set by the bath's temperature. 

To the lowest order in j* the state reduces to the 
near-equilibrium form pg^e ^"f^'^ +2n-(i-„») -^i with 
Peff = j[ In "'^^r and Ji — iQi the current operator 
with expectation (Ji) = j*. At high temperatures /3eff 
reduces to /3 = (/3l + Pr)/2. 

From this analysis, it appears that the identification 
(3eff — ao/h is physically grounded. For a finite current 
value one may use the symmetry property ai{—j*) — 
(— 1)'q!/(j*) to split J2i^o'^iQi ^ current-like part 
y — J2i>o '^2i+iQ2i+i , which is odd under boundary re- 
flection, and a remaining even part K, — X]/>i OL21Q21 [IBl- 
The steady-state is written then p*g ~ g-PesjH''-y-K. 
with an effective inverse temperature /3e// — ao/h that 
can be decomposed into /3e// — Pconf + A(j*^) where 
Pconf = In is the configurational (level popula- 
tion) part and A the current contribution. Since A ~ j*^ 
at small currents, the current contribution to the temper- 
ature does not show up in the linear regime. 

In summary, we have obtained the exact QNESS of a 
finite XX chain in contact at both ends with repeated- 
interaction reservoirs. We have shown that in the steady 
state, the system is completely specified by two quanti- 
ties, namely the magnetization profile (particle density) 
and the associated current. The flatness of the magneti- 
zation profile is related to the integrability of the model, 
leading to the violation of Fourier law since the system 
shows ideal conductivity. The (^NESS is given by the 
generalized Gibbs state e~^''ff^ -y-K. inverse tem- 
perature Peff with respect to a reference Hamiltonian 
. The many-body terms K, and y, build on system 
conserved quantities, are respectively symmetric and an- 
tisymmetric with respect to f . 
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